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Summary 


The  report  gives  an  account  of  a  mathematical  model 
for  estimation  of  the  ground  coverage  as  a  result  of  an 
attack  with  war  gas.  The  purpose  of  the  model  is  to  furnish 
an  estimate  of  the  size  and  distribution  of  the  ground  cov¬ 
erage  following  a  given  hostile  effort. 

Order  No:  101301 

Collaborators:  B.  Rydgren  (Appendix  1),  A.  Oehman  (Appen¬ 
dix  2),  A.  Linderot. 


Introduction 


The  purpose  of  this  work  is  to  study  a  model  for 
estimation  of  the  ground  coverage  after  an  attack  with 
liquid  war  gas.  For  estimation  of  the  extent  and  intensity 
of  thi  ground  coverage,  knowledge  of  the  existing  wind  and 
temperature  conditions  in  the  strata  in  question  is  of 
great  importance.  Works  concerning  the  influence  of  the 
meteorological  factors  upon  the  dispersion  are  described  in 
bibliographies  by  Sutton  *■)  and  Pasquill  2),  in  Sweden  such 
problems  have  been  dealt  with  by  Froessling  ^)»  4),  5), 
Aurivillius  Wedin  Wickerts  and  others.  In  refer¬ 
ence  7)  is  shown  a  model  for  estimation  of  the  amount 
required  for  an  attack  with-war  gas  in  aerosol  form,  which 
has  been  modified  here  to  make  it  applicable  to  gas  in 
liquid  form.  Reference  7)  also  gives  certain  guidelines 
for  estimating  the  ground  coverage  obtained  during  spraying 
of  liquid  war  gas,  but  the  influence  of  atmospheric  dis¬ 
turbances  upon  the  drops  during  the  spreading  is  not  taken 
into  consideration.  In  the  model  proposed  here,  this  factor 
is  taken  into  account. 


1.  General  Guidelines 


There  are  two  main  methods  of  spreading  war  gases  in 
liquid  form:  a  bomb  containing  the  gas  can  be  detonated  by 
means  of  a  relatively  small  explosive  charge,  or  the  gas  can 
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be  sprayed  from  tanks  attached  to  aircraft.  The  liquid  will 
be  in  the  form  of  drops  of  various  sizes,  most  of  which  will 
have  diameters  of  between  50  ^m  and  l,000^tm.  Even  though 
the  amount  by  weight  present  in  the  form  of  drops  less  than 
50  ft. m  in  diameter  is  insignificant,  this  part  of  the  sub¬ 
stance  is  of  great  importance  to  poisoning  of  the  atmosphere 
for  a  long  time  after  the  spreading.  In  general,  the  height 
from  which  the  spreading  is  done  varies  between  15  and  100 
meters.  During  the  fall  of  the  drops  toward  the  ground,  a 
vertical  speed  equilibrium  (vz),  depending  upon  drop  size 
and  the  density  of  the  liquid,  is  rapidly  reached.  The 
speed  of  the  drops  in  the  wind  direction  may  be  regarded  as 
equal  to  the  wind  speed  (vx) .  Together  with  the  original 
position  of  the  drop,  vz  and  vx  determine  the  expected  path 
of  the  drops.  Atmospheric  turbulence  will  cause  departure 
from  the  expected  path.  Here,  we  suppose  that  the  drops 
during  their  descent  spread  as  if  they  were  gas  particles 
(see  Section  53.).  The  departures  of  the  particles  can  be 
described  by  the  dispersion  values  or  <r  values  obtained  from 
meteorological  data.  These  depend  upon  existing  wind  condi¬ 
tions,  the  descent  period  of  the  drops,  and  upon  the  type  of 
atmospheric  stratificiation.  The  cr  values  in  the  wind  direc¬ 
tion  (<rx),  at  right  angle  to  the  wind  direction  (<r  y),  and 
vertically  (CTz)  are  here  assumed  to  be  equal  (see  Appendix 
1).  This  is  a  rather  exact  approximation  for  indifferent 
stratification,  and  often  for  stable  stratification.  It  is 
at  these  types  of  stratification  that  a  scattering  of  this 
kind  can  be  expected  to  occur,  while  the  effect  at  a  labile 
stratification  might  be  very  hard  to  judge.  The  scattering 
device  may  also  locally  increase  the  atmospheric  turbulence. 


2.  The  Primary  Cloud 

The  primary  cloud  is  the  designation  for  that  cloud 
of  liquid  drops  which  forms  after  the  spraying.  It  presents 
various  appearances,  depending  upon  the  production  method. 

If  the  primary  cloud  is  produced  by  means  of  a  bomb,  it  is 
often  of  spherical  or  similar  shape;  if  formed  by  means  of 
spraying  from  a  tank  it  appears  as  a  long  string  of  rela¬ 
tively  little  height  and  width.  The  primary  cloud  is  assumed 
to  be  built  up  of  a  number  of  right-angle  parallelepipeds, 
with  sides  1*,  ly,  and  lz.  These  elements  are  called  cloud 
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elements.  The  distribution  of  drop  sizes  within  the  various 
cloud  elements  is  in  the  following  assumed  to  be  generally 
the  same  in  all  elements,  dcch  cloud  clement  is  thought  to 
be  built  up  of  a  number  of  stratified  unit  elements,  each  of 
which  contains  a  narrow  size  fraction  of  the  drops. 

For  estimation  of  the  coverage  from  a  primary  cloud, 
the  cloud  is  assumed  to  be  divided  into  unit  elements.  The 
coverage  from  each  unit  element  is  calculated  and  adding 
these  coverages  gives  the  coverage  from  the  primary  cloud 
(see  61).  Programming  for  a  /data  processing.7  machine  has 
been  worked  out  for  such  calculations;  it  is  shown  in  Appen¬ 
dix  2.  In  this  case,  then,  the  three  dimensions  of  the 
primary  cloud  and  the  drop  size  are  nddea  up.  In  certain 
cases  (see  62)  integrals  can  be  substituted  for  these  sums, 
which  lessens  the  time  required  for  the  calculations. 


3.  Distribution  of  Drop  Sizes 

Drops  of  liquid  of  various  sizes  are  produced  by  the 
formation  of  the  primary  cloud.  The  distribution  of  the  drop 
sizes  depends  upon  the  method  of  delivery.  Among  the  dis¬ 
tribution  functions  available  in  the  literature,  the  so- 
called  Nukiyama-Tanasawa  distribution  has  been  selected 
here: 


dN  vP 
—  ■  a.XK.exp 

ax 


(3.1) 


In  practical  cases,  it  may  become  necessary  to 
examine  pragmatically  the  validity  of  this  assumption.  In 
particular,  it  may  be  important  to  check  the  distribution 
of  very  small  and  very  large  drops,  in  the  former  case 
because  of  the  increased  risk  of  poisoning  via  the  respira¬ 
tory  system,  and  in  the  latter  case  because  of  the  greater 
mass. 


In  (3.1)  a,  b,  p,  and  q  are  parameters.  N  repre¬ 
sents  the  relative  number  (fraction  of  the  total),  and  X  the 
drop  diameter  in  fixa,  Thus  the  fraction  of  drops  to  which 
XX<X  «£  X2  apply  is  x 

/dN 

X, 


If  parameters  p  and  q  are  assigned  values  2  and  1, 
respectively,  good  agreement  is  found  with  the  empirically 
found  distribution  curves  of  reference  5). 

The  volume  of  all  drops  having  a  diamter  of  4  X 
becomes: 

V  ./*  JL  .  X3  .  dN  (3.2) 

o  6 

From  (3.1)  and  (3.2)  and  with  the  figure  values  for  p  and  q 
there  is  obtained: 

dV  „  a  .  %  .  X5  .  exp  {-bx}  (3.3) 

dX  6 

The  following  designations  are  introduced: 


(3.4) 


in  which  VQ  designates  the  total  spray  volume. 

The  various  central  values  in  question  for  this  distribution 
are  obtained  as  follows: 


(3.5) 


gives  X«  D,  where  D  designates  the  maximum  value  or 
the  type  value  (volume-mode -diameter)  for  the  frequency 
function  dV 
dX 
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The  parameters  a  and  b  can  now  be  expressed  In  V0  and  D  as 
follows: 


a  ■ 


6.b6.V( 
*  .5! 


(3.7) 


(3.8) 


Further,  there  Is  the  following  connection  between  D  and  Dm 
when  p  ■  2  and  q  •  1: 

Dm  •  D. 1.1340 

which  is  obtained  from  (3.3),  (3.6),  (3.7)  and  (3.8). 

If  (3.7)  and (3.8)  are  introduced  into  (3.3)  there  is 
obtained 


dV 

dX 


55.V 


24.  D* 


•exp  {■  x' x} 


(3.9) 


In  direct  spraying  from  aircraft  there  is  obtained,  by  means 
of  a  formula  according  to  Merrington  and  Richardson 

h' 


Dm  in  Km  from  the  relation 


0.57.107.  Vv 


(3.10) 


a 


in  which  V  ■  kinematic  viscosity  in  Stokes 


va  •  relative  speed  between  the  liquid  jet  and 
and  the  ambient  air,  cm/sec. 

For  devices  tested  in  aircraft  experiments  by  Froessling 
and  Aurivilius  5)  there  was  obtained 


io7  .  V7 
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Din  a* 


(3.11) 


In  the  division  of  a  cloud  element  into  unit  elements, 

(3.9)  is  approximated  with  a  •■•to  frequency  function. 

The  X  axis  is  divided  tine  -n  nuul-o.  of  intervals  of  equal 
size.  The  intervals  become  (i-l).b’  _o  i  b’  where  i  assumes 
the  values  1,  2  ...  ta.  In  the  discrete  frequency  func¬ 
tion,  X  can  only  assume  the  value  corresponding  to  the  m 
number  of  the  middle  points  of  the  intervals,  thus 
Xi  *  . 5.b*  ,  X2  -  1.5  .  b‘  ,  Xk  =  Oc  -  .5).  b»  .  .  .  X,*  « 

(m  -  .5).  b* .  In  the  interval  around  X<K  the  spread  volume 
becomes 

k  .  b’ 

J  dV 

(k  1)  .  b’ 

m  is  selected,  for  example,  so  that 
m  b 

f  dV  >  0.959.Vq 
o 


4.  The  Speed  of  Descent  of  the  Drop" 

The  speed  of  descent  (vz)  of  the  drops  depends  on 
their  diameters  and  density.  Here,  as  in  reference  7), 
the  following  relatinn  between  drop  diameter  ( /nn)  and 
speed  of  descent  (m/s)  is  used: 

/Ci  ♦  C2  los  X  +  c3  (logX)£7 

Vz  •  k.  X  ~  (4.1) 

in  which  k  «  Stokes'  constant  -  3.0044  .  10"^  in  air 
at  70°F  (21.1°C)  for  a  drop  density  of  1  g/cm-*. 

For  drops  of  X  <  53  ^m,  Ci  •  2  and  C2  •  0. 

For  drops  in  the  interval  53  <  X  <  800  pm,  * 

2.01302,  C2  •  .12055,  and  =  -  .07422 

Since  the  number  of  drops  of  diameter  7  800  /nn  are 
probabl  very  rare,  (4.1),  with  the  latter  values  for  the 

-  6  - 


constants,  has  been  used  in  the  area  X  800  m.  nxact 
values  may  be  obtained,  for  example,  according  to  Froess- 
ling  10>. 

It  is  assumed  here  that  evaporation  from  the  drops 
is  negligible.  When  this  is  not  the  case,  however,  the 
calculations  can  be  based  on  the  actual  time  of  descent. 


5,  The  Spread  of  the  Unit  element  and  Coverage  From  It 

51.  Designations  (see  Fig  1) 

x,  6,  and  z  designate  the  positions  in  a  right 
angle  coordinate  system  with  the  x-axis  in  the  wind  direc¬ 
tion.  The  positive  z-axis  designates  height  above  the  ground. 
The  y-axis  is  at  a  right  angle  to  the  x-z  plane. 

(Xi  ;  y^)  position  in  the  ground  plane. 

t  variable  time.  t*0  is  time  of  formation  of  the  cloud 

t8  point  of  time  when  protection  is  obtained 

r  distance  from  a  point  to  the  center  of  the  cloud 

C  concentration 

Cm  concentration  at  the  center  of  the  cloud 

(xnj  ;  ym  ;  z^j)  position  of  the  center  of  the  cloud 

(xo  J  Fo  I  z0>  positions  of  the  center  points  at  t  -  0 

rQ  distance  from  the  center  of  the  c loud  to  a  point  where 
the  concentration  is  .05  .  Cm 

t  i  point  in  time  when  the  distance  from  a  point  (x^  ,  y^) 
to  the  center  of  the  cloud  is  minimal 

r0^  is  r0  at  point  in  time  tQ^ 
vx  median  wind  speed  (see  Appendix  1) 
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vz  speed  of  descent 

w  the  amount  of  matter  in  the  cloud 

Ax*  Ay*  Az  constants  for  the  turbulent  diffisulon  in 

x-,  6-  and  z-  links 

<r„  <r„  <r„  the  spread  of  drop  around  its  expected  path 

x>  y>  z  in  X-;  y-  and  a-  links 

<r  the  spread  of  a  drop  about  the  expected  path  of  the 
center  of  the  cloud 

n  layering  constant 

D(xi  ;  dosage  at  the  point  (x-^  ;  y^)  up  to  the  time  t8 
B(xi  ;  y^)  coverage  at  the  point  (x^  ;  y^)  up  to  the  time  ts 


52.  The  Development  of  the  Unit  Elements 

At  the  point  of  time  t  *  0  the  unit  element  has  the 
sides  lx,  ly  and  1  and  contains  drops  of  so  narrow  a  size 
fraction  that  the  drop  size  may  be  regarded  as  constant. 

The  concentration  at  all  points  of  the  cloud  is  assumed  to 
be  the  same  at  the  point  of  time  t  =  0.  Thus,  if  one  regards 
the  concentration  distributions  in  x,  y,  or  z  links  there  is 
obtained  a  rectangular  concentration  distribution  with  the 
variation  widths  lx,  ly.9l„  respectively.  The  variance  for 
such  a  distribution  is7ivZ.  Immediately  after  the  point 
of  time  t  *  0  the  particles  become  exposed  to  the  disturb¬ 
ance  caused  by  atmospheric  turbulence.  Due  to  this  dis¬ 
turbance,  each  drop  deviates  randomly  from  its  expected 
path.  The  deviations  from  the  expected  path  in  both  x,  y, 
and  z  links  are  regarded  as  normally  distributed.  Thus, 
the  distribution  of  the  particles  around  the  path  of  the 
center  of  gravity  of  the  cloud  will  be  composed  of  a  rec¬ 
tangular  distribution  (the  original  position  of  the  drop  in 
the  unit  element)  and  a  normal  distribution  (atmospheric 
disturbance).  The  spread  of  the  normal  distribution,  which 
grows  with  time,  after  a  rat.ier  short  time  becomes  large 
in  relation  to  the  constant  spread  of  the  rectangular  dis¬ 
tribution,  The  concentration  of  the  cloud  can  then  be 
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wind  direction. 


regarded  as  determined  by  a  normal  distribution,  in  which 
the  spread  is  composed  of  tha  spread  oJ  the  two  distribution. 

As  previously  mentioned,  the  expected  path  of  the  cen¬ 
ter  of  gravity  of  the  unit  element  is  determined  by  the 
original  position  of  the  center  of  gravity  and  wind  speed 
and  speed  of  descent.  The  wind  speed  usually  increases  with 
the  altitude  (the  wind  speed  profile),  and  the  path  there¬ 
fore  is  not  a  straight  line.  At  the  altitudes  we  are  con¬ 
cerned  with  here  the  spread  values  have  been  assumed  to  be 
Independent  of  the  altitude.  Further,  it  has  been  assumed 
that  a  designed  so-called  equivalent  wind  speed,  vx,  a  con¬ 
stant  in  the  strata  in  question,  can  he  substituted  for  the 
wind  speed  increase  with  altitude  (see  further  Appendix  1). 
This  wind  speed  equivalent,  vx,  depends  for  a  certain  ground 
wind  speed  and  layer  constant  upon  the  original  altitude  of 
the  unit  element.  The  wind  equivalent  is  so  selected  that 
the  expected  point  of  impact  of  the  center  of  gravity  is  the 
same  as  if  the  real  wind  speed  profile  had  been  used.  In 
accordance  with  this  simplification  the  path  of  the  center 
of  gravity  becomes  a  straight  line.  Thus  the  center  of 
gravity  is  assumed  to  be  displaced  at  the  constant  speed 
vx  in  the  wind  direction,  while  at  the  same  time  it  falls 
at  the  constant  speed  vz.  The  displacement  x^  -  xff  in  the 
wind  direction  during  the  fall  from  the  original  altitude 
z0  becomes 


x 


1  " 


z 


o 


V 


z 


53.  Calculation  of  the  Ground  Coverage  From  a  Unit  Element 

The  unit  element  is  considered  to  become  a  cloud  with 
the  concentration  gradient  determined  by  the  normal  distribu¬ 
tion  long  before  the  ground  level  is  affected  by  the  cloud. 

It  can,  therefore,  be  assumed  that  the  unit  element  immediate¬ 
ly  turns  into  a  normal  distribution  cloud  with  the  variances 
12X/12,  l2  /12,  and  12Z/12  in  the  x-,  y-,  and  z-links, 
respectively.  After  the  point  of  time  t  ■  0  the  spread  of 
the  cloud  caused  by  atmospheric  disturbances  begins.  The 
variance  of  gas  element  around  its  expected  path  due  to 
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this  cause  is,  according  to  Sutton  (see  Appendix  1); 


or  2 

X 

=  0.5  • 

A2  . 

X 

(t 

c 

I 

CM 
/— s 

.>* 

• 

CM 

>* 

b 

=  0.5  • 

A2y  • 

(t 

-  v  2-n 
*  vx> 

-Z2 

=  0.5  • 

A2z  • 

(t 

•  vx>2"n 

(5.1) 


For  altitudes  of  <100  m  it  is  often  approximately 
holds  true  for  stable  and  indifferent  stratification  that: 

Ax  =Ay»Az  =  A 
and  therefore 


It  is  here  assumed  that  the  drops  during  their 
descent  have  the  same  deviations  ns  gas  elements  from' their 
expected  paths,  so  that  the  <r  values  according  to  (5.1)  are 
valid  regardless  of  the  drop  size  (W.  Schmidt  1:L5 .  For  large 
drops,  however,  diminished  spread  may  occur  for  two  reasons: 
the  restraining  effect  of  inertia,  and  because  the  drops 
enter  new  urbulence  more  quickly  because  of  their  greater 
speed,  (see  Pasquill  ',  pp  112,  226).  In  certain  cases  of 
calculations  according  to  the  formula  below  it  may  be  neces¬ 
sary  to  employ  a  values,  which  in  this  respect  vary  with 
drop  sizes. 

If  one  selects  lx  •  1  =  ly  =  1  the  variance  of 

the  drops  around  the  expected  path  of  the  center  of  gravity 
for  the  three  main  direction;;  becomes: 

a2  •  0.5  A2  •  (t  •  vx)2"n  12/12  (5.2) 

Since  the  variances  in  the  main  directions  are 
assumed  to  be  equal,  points  of  the  same  concentration. are 
located  on  a  spherical  shell.  At  tnc  distance  r  from  the 
middle  point  of  the  cloud  this  concentration  is 

C(r)  -  •  exp  {-r2/2<r2} 
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Here  C(rQ)  -  0.05  •  Cm 


°-05  cm  '  cm  •  exP  K  2/?'2} 

which  gives 

ro2-6.0-<r2  (5.3) 

Thus 


C(r)  »  Cm  .  exp  |-3r2/r02] 


(5.4) 


C_  can  be  determined  from  the  amount  of  substance  in  the  unit 
element . 


w  =  4  •  J  r  •  C(r)  •  dr 


which  gives 


w  .  3  .  /3 

31  .  ft. 

r  O 


(5.5) 


If  (5.5)  is  introduced  into  (5.4)  there  is  obtained 

C(r)  =  ”  '  3  '  £L  exp  j-3r2/r02  j  (5.6) 

X  .  ■{%  .  r3  * 

'  o 

The  position  of  the  middle  point  at  time  t  is  deter¬ 


mined  by  the  following  equations: 

*m  =  xo  +  c  *  (5*7a> 

ym  *  y0  (5.7b) 

zm  "  zo  *  t  •  vz  (5.7c) 

The  distance  from  the  points  (x^  ;  y^)  to  the  middle 
point  of  the  cloud  can  be  witter: 

r2  *  <X1  ”  Xm^2  <yi  '  ym)2  (5-8> 


12 


If  equations  (5.7)  are  introduced  into  (5.8),  r2  is 
obtained  as  a  function  of  the  time 

r2  *  (xx  -  x0  -  t  •  vx)2-t-  (yi  -  y0)2+  (zo  "  c  *  vz)2  (5‘9> 

If  (5,9)  is  introduced  in  (5.6)  the  concentration  in 
points  (x^  ;  y]_)  is  obtained  as  a  function  of  the  time. 


(yi  -  y0>2  +  <zo 


t  •  vx)2  + 


(5.10) 


For  calculation  of  the  dosage  ano  the  coverage  a 
simplification  is  introduced:  the  cloud  is  assumed  to  be 
of  constant  size  during  passage  of  the  points  (x^  ;  y{). 
During  this  passage,  the  cloud  is  assumed  to  have  the  Cra 
and  r0  values  which  exist  when  the  distance  from  the  points 
to  middle  point  of  the  cloud  is  minimal.  (Before  and  after 
this  point  the  cloud  is  in  reality  smaller  and  larger, 
respectively.)  The  minimal  distance  occurs  at  time  t 
The  radius  of  the  cloud  is  then  r  ,  ■  rQ  (tQ^)  •  t 
is  obtained  through  reducing  (5.9),  which  results  in 


vx(*l  -xQ)  - 


(5.11) 


Introduction  of  t  ^  in 


(5.2)  and  (5.3)  gives: 

2-n  9 

<*olV  +  °-5  '  1 


If  1  m  is  selected  for  the  size  for  1,  the  last  term  can  be 
eliminated.  The  dosage  in  points  (x^  ;  y^)  is  obtained  from: 


D(xx  ;  yx) 


r 


j 

o 


C  ( t )  •  c.  t 
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If  (5,10)  is  Introduced,  there  Is  obtained: 


!>(*!  ;  yx)  * 


1  rol 


w  .  3  ./3  exp  (-  — ll  (yi  -  y0)zJ 


*  •  if*  .  rol3 


I  exp  {'  “^7  •  [<*1  -*0  *  t-vx)2+  (t-'i  -  *0>2ji  *  <: 

o  rol  ' 

Through  the  substitutions 

E  -  vx  (xt  -  x0)  +  vz  .  z0 

A2  =  vx2+  Tl* 

and  the  transformation 

~t  (*t  -  V  *  “2 

rol2.  *  2 


there  is  obtained 
w  •  3 

D(*l  ;  yi)  - 


«*P  {  -  -i—  [(y^  -jjK  (xx  -  *0)2+z, 
*  rol 


I  Jr 


-  »|2 


rol  • 


A  .  * 


®  2  .  du 


(5.13) 


where  the  integration  limits  b  and  a  are  defined  by: 


b  «  (A*t  -  f) 
rol  A 


a  *  - 


/6  •  E 


rol  *  A 


.12) 
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The  coverage  (quantity/m2)  is  in  simple  relation  to 
the  dosage  (time  •  quantity/m-^) ,  which  is  apparent  from 
the  following:  The  additional  quantity  dB  during  the  time 
interval  t  for  a  point  at  ground  level  derives  from  drops 
found  between  ground  level  and  the  height  At  •  vo  above 
ground  level.  The  concentration  during  the  time  interval 
d  t  in  the  area  of  the  point  is  constant  and  equals  C(t), 

The  additional  quantity  can,  therefore,  be  expressed: 
dB  »  At  •  vz  •  C(t) 

After  Limes  transition  (to  differentials)  and  inte¬ 
gration  there  is  obtained: 

tg 

B(xi  ;  y^)  -  vz  •  J  C(t)  dt  =  vz  •  D(x^  ;  y^) 


The  coverage  can  thus  be  obtained  from  the  dosage 
by  multiplication  by  the  speed  of  descent. 


‘ol 


Summary  of  Formulas 

Vx  (*1  -  *0)  +  vz  .  z0 

2  2 
vx  ♦  vz 

VX  (X1  *  *o>  vz  *  ze 


-2 


E 


D(xt 


3  •  A2  (tol  •  vx)2-"  *  0.5  .  I2 

3-wexp  -|_i_  J<yi  -  ya)2  ♦  (*j  -  x0)K  . 

*i>  8 _ roi _ _ • 

x  .  A  .  rol2 
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«  =  .  V*  •  8 

rol  •  A 

b  •  j£  (A  t,  -  I) 
rol  A 

B  (xt  ;  yj)  ■  v2  .  D(x!  ;  yj) 

6.  Calculation  of  the  Coverage  From  a  Primary  Cloud 
61.  Addition  of  Unit  Elements. 

The  primary  cloud  is  described  by  the  positions  and 
contents  of  the  cloud  elements.  The  middle  point  of  the 
i-te  cloud  element  is  designated  (xQ.;  y  r  .).  The  total 
amount  of  substance  in  the  primary  ctoudis  designated  W. 

The  amount  of  substance  in  unit  element  j  of  the  1  cloud 
element  is  designated  .  Therefore: 

j  vu  *  vi 

if  wu  "i'1  ’  "  • 

in  which  w^  designates  the  amount  of  substance  in  the  1-te 
cloud  element. 

The  drop  sire  distribution  within  the  various  cloud 
elements  is  assumed  to  be  the  same.  In  the  preceding  section 
the  coverage  from  a  unit  element  was  calculated.  The  cover¬ 
age  from  a  primary  cloud  is  obtained  through  addition  of  the 
coverages  from  all  of  the  unit  elements.  These  calculations 
are  carried  out  by  means  of  computer.  A  program  for  the 
computer  is  shown  In  Appendix  2. 
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In  dealing  with  a  primary  cloud  originating  In  spray¬ 
ing  from  a  tank  attached  to  an  aircraft,  the  primary  cloud 
can  with  adequate  exactness  be  described  as  a  number  of 
cloud  elements  lying  above  one  another  and  with  their  cen¬ 
tral  points  lying  in  a  straight  line  in  space.  The  size, 
form,  and  position  of  the  primary  cloud  is  thus  determined 
by  the  position  of  the  first  cloud  element  (x0j  ;  y0i;  *0i>, 
the  number  of  cloud  elements,  M,  and  the  position  of  the  last 
cloud  element  (x0^;  y^;  Zqm).  The  position  of  1  cloud 

element  is  then  (x0£j  y0£j  zQ\)  in  which 

*01  *  *ol  *  (I'D  — 

M-i 

and  the  corresponding  expressions  for  y0£  and  zQ^. 

Figures  2  shews  the  ground  coverage  calculated  for  a 
unit  element,  produced  at  an  altitude  of  25  m,  having  sides  of 
1  m,  with  wind  speed  at  7  tn/s,  and  the  layering  constant,  n, 
at  .35,  i.e.,  A -.07  (see  table  in  Appendix  1).  The  substance 
content  of  the  unit  element  is  1  g  and  the  substance  density  1. 

Figure  3  shows  under  the  same  conditions  B(x;  0),  i.e. 
the  coverage  on  the  x  axis,  from  unit  elements  corresponding 
to  two  different  drop  sizes. 

Figure  A  shows  the  ground  coverage  from  a  primary 
cloud  sprayed  at  right  angles  to  the  wind  direction  at  the 
altitude  of  25  m.  The  cloud  has  a  length  of  1,000  m  and  is 
regarded  «s  being  built  up  of  a  row  of  cloud  elements,  each 
of  which  has  sides  of  1  ro.  Wind  speed,  layering  constant, 
and  substance  density  as  in  Figure  1.  The  cloud  has  a  drop 
size  distribution  according  to  (3.9),  with  the  median  dia¬ 
meter  at  280  The  substance  content  is  300  g/tn  ^sic; 

should  probably  read  300  g/mJ7.  In  the  calculation  of  the 
coverage  the  cloud  elements  were  divided  into  17  unit  ele¬ 
ments,  each  of  them  including  a  drop  size  interval  of  50  /‘m. 

An  investigation  of  the  effect  of  interval  length  upon  the 
coverage  has  shown  that  a  finer  division  than  50  /*m  is  un¬ 
necessary.  The  substance  content  and  the  speed  of  descent 
for  the  various  unit  elements  are  shown  in  the  table  below. 

The  speed  of  descent  of  a  unit  element  is  calculated  from 
(4.1),  in  which  the  drop  size,  X,  is  equal  to  the  middle 
point  of  the  Interval. 


*> 


Drop  size  for  the 
unit  cube  (pm) 

Substance 
content,  w(g) 

Speed  of  descent 
v*  (m/s)  • 

0-50 

0.190 

0.019 

50  -  100 

5.055 

0.154 

100  -  150 

21.079 

0.350 

150  -  200 

40,496 

0.578 

200  -  250 

51.682 

0.820 

250  -  300 

51.386 

1.067 

300  -  350 

43.211 

1.314 

350  -  400 

32.228 

1.557 

400  -  450 

21.966 

1.794 

450  -  500 

13.958 

2.025 

500  -  550 

8.385 

2.247 

550  -  600 

4.812 

2.462 

600  -  650 

2.658 

2.668 

650  -  700 

• 

1.421 

2.866 

700  -  750 

0.739 

3.055 

750  -  800 

0.375 

3.237 

800  -  850 

0.360 

3.411 

In  chi*  case,  where  Che  dlssemlnaClon  was  done  aC 
rlghc  angle  co  Che  wind  dlrecclon  and  aC  a  constant  altitude, 
lc  was  possible  Co  lncroduce  certain  simplifications  into 
Che  computer  program.  As  indicated  In  the  foregoing  and  as 
Is  apparent  from  Appendix  2,  one  can  by  mean*  of  the  computer 
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program  worked  out  obtain  for  an  arbitrary  point  on  the 
ground  level,  the  coverage  from  a  primary  cloud  element 
consisting  of  a  cloud  element  in  an  arbitrary  configuration 
in  space.  By  means  of  a  minor  adjustment  to  the  computer 
program,  it  is  in  addition  possible  to  vary  the  drop  size 
distribution  from  cloud  element  to  cloud  element. 


62.  Primary  Cloud  Integration 

In  61  the  coverage  at  a  point  (x^  ;  y^)  was  calcu¬ 
lated  by  summation  of  (5.13)  on  drop  size  X  and  on  the 
dimensions  x«,  yQ  and  zQ  of  the  primary  cloud.  The  total 
coverage  at  (xj  ;  y^)  can  thus  be  written 


\  \  \  £  vz  •  D<*1  ;  Yi) 

o  yQ  z0 

If  the  form  of  the  primary  cloud  can  be  given  in  a 
simple  mathematical  expression,  or  more  of  the  symbols 
of  the  above  formula  can  to  advantage  be  replaced  by  inte¬ 
grals.  Froessling  has  shown  such  calculations  for  gas 
clouds. 


Here  we  shall  only  show  that  the  primary  cloud  can 
be  regarded  as  a  linear  cloud  at  altitude  z0,  where  the  angle 
between  the  direction  of  spraying  and  the  wind  direction  is 
determined  by  the  constant  k  in  relation  xQ  :  k  •  yQ. 

The  coverage  at  (x|  ;  y^)  is  now  calculated  from: 

B(xx  ;  yi)  •  L  v  J  D(xt  ;  yi)  dy 

X  v 

'  o 

yo  ;  y^)  •  dyQ  is  the  dosage  originating  from 

yQ  drops  of  a  certain  size,  giving  the  speed  of  descent  vz. 
This  integral  will  probably  be  of  interest  in  many  appli¬ 
cations,  and  the  solution  will  therefore  be  shown. 
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20 


250  « 


9 


f  d(*i  ;  yi)  dy0 
y0 


where  wy  is  the  amount  of  substance  per  unit  of  length 
in  the  y  link,  w  is  calculated  from  the  amount  of  substance 
per  unit  of  length  in  the  cloud,  w^,  with  the  relation 


0L-  k2  +  1 


xx  •  k  +  yx  - 


k  vx  Gi  •  X1  "  vz  •  z0) 


1  :  JL.  .<3e.fi  -j.) 


rol 


r0^  is  taken  from  the  formulas  for  rQn  and  t_^,  page  15, 
with  xQ  •  k«y^.  This  relation  is  used  also  in  calculating 
the  integration  limits  a  and  b. 

From  the  last  formula,  the  integration  limits  and 
q2  are  calculated  from  the  corresponding  limits  for  y0. 
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In  formula  (6.1)  it  is  possible  to  remove  the  paren¬ 
theses  from  the  integral  and  to  give  it  a  median  value. 
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APPENDIX  1 

THE  EFFECT  OF  METEOROLOGICAL  FACTORS  UPON  THE  MODEL 
Bo  Rydgren,  Staff  Meteorologist 

To  make  the  calculation  method  described  in  the  fore* 
going  applicable  in  the  simplest  way,  the  wind  speed  vx, 
variable  with  altitude,  has  been  replaced  by  a  so-called 
"equivalent  wind,"  vx,  which  gives  the  middle  points  of  the 
individual  drops  the  same  point  of  impact  as  if  vx  were 
used, 

vx  is  obtained  from 


/  V 


.  dt 


(Bl.l) 


in  which  t^  is  the  time  from  the  moment  of  release 
to  the  time  the  particles  reach  the  ground. 

According  to  Deacon  12)  the  real  wind’s  dependence  on 
the  altitude  can  be  written 


vx(z)  =  k 


(B1.2) 


in  which  z  is  the  altitude  at  which  the  wind  speed  is 
vx,  n  the  layering  constant  (see  below),  and  k  a  constant, 
the  value  of  which  can  be  determined  from  a  wind  speed  mea¬ 
sured  at  a  certain  altitude.  If  for  example,  the  wind  at 


the  altitude  of  5  m,  v(5)  is  assumed  known,  vx  is  obtained 
from  (B1.2)  as  follows 


log  v„(z)  =  -0-  (log  z  -  log  5)  ■+■  log  vx  (5) 
2-n 


(B1.3) 


As  is  seen,  wind  direction  variation  with  altitude 
has  not  been  dealt  with,  although  in  certain  cases  it  may 
have  a  considerable  effect  upon  the  spreading  (due  to  con¬ 
vection,  ground  friction,  or  in  connection  with  speed  of 
descent  of  the  drops,  for  example). 

If  the  speed  of  descent  of  the  drops  is  assumed  to  be 
constant,  (Bl,l)  can  be  written: 


vx  *  2o 


•  / 


vx  •  dz 


(B1.4) 


in  which  zQ  is  the  release  altitude. 

From  (B1.2)  and  (Bl.4)  there  is  obtained: 

2 


x  z. 


k.  .  z0 


(B1.5) 


Finally,  from  (B1.2)  and  (Bl.5)  there  is  obtained  for  the 
altitude  Zw  with  a  calculated  k  value  for  the  measured  wind 
speed  v(z^; :  n 


v*  -  a  -  ->  •  vx(*h)  •  ( 


(B1.6) 


From  (Bl.6),  vx  can  thus  be  calculated  for  a  given  layering, 
represented  by  the  constant  n  (see  below)  and  for  known  wind 
speed  v  (z.  ),  measured  at  the  altitude  z^  (see  equation 
B1.3).  *  " 
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It  should  be  added  that  the  above  reasoning  can  be 
applied  only  for  altitudes  up  to  about  100  m,  and  that  the 
wind  speeds  calculated  from  (B1.2)  should  be  regarded  only 
as  probable  values,  from  which  great  deviations  may  occur  in 
certain  cases. 

According  to  Sutton^),  the  following  connections  can 
be  used  for  determination  of  the  time  dependency  for  the 
variance  in  the  particle  distribution  within  the  individual 
drop  cloud: 

"x2  *  °-5  42x  2'n 

%2  -  0.5  ■  t)  2‘"  (B1.7) 

2  n  c  2  x—  .2-n 

(T  *  0.5  A  (v  .  t  ) 

2  Z  * 

where  the  diffision  constants  Ax,  A  ,  and  have 
been  determined  for  various  wind  speeds  in  field  experiments. 
In  this  case,  the  previously  defined  "equivalent  wind"  vx 
was  used. 

It  can  be  shown  that  at  indifferent,  and  often  at 
stable  layering,  within  0  -  100  m  the  following  is  valid 

A  x  Ay  ~  ^ 

For  calculation  of  the  successive  spread  of  the  drop¬ 
let  cloud  from  the  moment  of  release  until  contact  with  the 
ground  the  A  values  according  to  the  table  below,  calculated 
for  various  layering  conditions  and  wind  speeds  are 

U3ed, 


Table,  Numerical  A  values  for  various  values  of 
v  and  n,  applicable  to  slightly  rolling 
terrain  and  spread  tiroes  4  15  min. 


ITiTTW 

0) 

vindhastighet ,  t 

'x 

m 

i 

2 

3 

n 

5 

- 1 

6 

B 

8 

9 

10 

r 

labil 

<  0,20 

0,52 

0,33 

0,28 

C,23 

0,20 

0,18 

0,15 

X 

- 

- 

neutral 

0,20- 

0,30 

0,21 

0,19 

0,1? 

0,15 

0,13 

0,12 

0,11 

0,11 

0,12 

0,13 

atabil 

0,30- 

O.bo 

0,10 

0,09 

o 

o 

03 

0,07 

0,06 

0,07 

0,07 

0,07 

0,08 

0,08 

© 

■yckat 

•tabu” 

>  0,40 

0,05 

-*■ 

o 

«*■ 

o 

1 

1 

10,03 

L _ 

0,02 

*_  i 

1 

- 

x  Probably  constant  in  the  continuation,  with  the 
wind  within  the  interval  £10  m/s. 

**  Rvalues  for  very  stable  layering  are  uncertain. 

/Legend7:  (1)  Type  of  layering;  (2)  Layering 
constant;  (3)  Wind  speed  v  m/s;  (4)  Very 
stable. 


APPENDIX  2 

DATA  MACHINE  PROGRAM  FOR  CALCULATION  OF  GROUND 
COVERAGE  FOR  SPREAD  OF  SUBSTANCES  IN  LIQUID  FORM 

Engineer  Arne  Oehraan 

The  program  shown  below  was  prepared  at  FOA  Section 
322  at  the  request  of  FOA  1.  The  program  is  coded  in  FORTRAN 
IV  to  be  run  through  the  FOA  data  procesiing  machine  IBM  7090/ 
1401. 


Basis,  In  and  Out  Values 


The  area  of  operation  is  a  right  angle  coordinate 
system,  in  which  x-y  designates  the  ground  plane  with  the 
x  axis  in  the  wind  direction  and  z  the  height  above  the 
x-y  plane. 

At  the  moment  of  spraying  the  sprayed  substance  is 
in  the  form  of  a  primary  cloud  build  up  of  a  number  of  cloud 
elements  of  the  same  shape,  Each  cloud  element  is  in  turn 
built  of  superimposed  unit  elements,  with  unit  elements  of 
the  same  type  assumed  to  contain  droplets  of  constant  size. 

The  division  into  unit  elements  is  described  in  the 
chapter  SUBROUTINE  MOFHST. 

In  values: 


NRS 

Operation  No. 

SN 

Layering  type 

V5 

Wind  speed  at  ground  level 

m/s 

NAMAL 

Number  of  targets 

NEMT 

Number  of  cloud  elements 

EL 

Cloud  element  side  ^ 

m 

DST 

Drop  size  interval  in  the  unit  element 

H  m 

DMM 

Mass  median  diameter 

W 

Amount  of  sub stance /cloud  element 

g 

XM,  YM 

Target  coordinates 

m 

XO,  YO, 

ZO  Cloud  element  coordinates 

m 

TS 

Time  to  protection  from  target 

s 

1)  The  cloud  element  is  assumed  to  be  cube-formed. 

Out  Values: 

Tabic  of  division  Into  unit  elements 
Coverage  on  targets 

Calculation  Procedure: 

1.  Data  entered 

2.  Median  speed  VM  arc  calculated  for  ail  heights  ZO 
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(SN/(2-SN)) 


VM  -  A-~  SN-  .  V5  .^ZOj 

3.  DELTA  2  «  3  •  t?"  is  calculated  for  all  median 
wind  speeds  VM. 

4.  The  cloud  elements  are  divided  into  unit  elements. 

/  Number  of  unit  elements /cloud  elements  s  NF.  Speed  of 

descent,  VZ,  and  share  of  substance,  DW  for  unit  elements 
are  determined.  Substance  density  assumed  to  *  1. 

5.  For  each  target  ground  coverage  at  time  TS  is 
calculated  through  adding  together  coverages  from  the  NEMT 
cloud  elements.  The  contribution  from  a  cloud  element  is  in 
turn  obtained  by  adding  together  the  contributions  from  the 
NF  unit  elements  in  the  cloud  element.  The  total  average 
for  a  target  is: 

NEMT  NF 

B^tot  *  Z  Z  BLG  (J,K) 

J-l  K-l 

in  which  BLG  (J,I<)  is  the  contribution  from  the  K 
unit  element  in  the  J  cloud  element. 

The  rest  of  the  calculation  procedure  is  apparent  from 
the  flow  scheme. 


SUBROUTINE  MOFHST 


Purpose  of  the  Routine 


The  routine  is  employed  for  dividing  a  cloud  element 
into  unit  elements  and  for  calculating  the  speed  of  descent 
of  the  respective  unit  elements. 


The  division  is  made  into  intervals,  determined  by 
the  parameter  DST,  and  the  speed  of  descent  is  calculated  for 
the  median  value  of  the  droplet  size  within  the  unit  element. 
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Basis.  In  and  Out  Values 


Call:  CALL  MOFHST  (DST,  DMM,  I,  DEL,  VZ) 

In  values:  Mass  median  diameter  (DMM)  ft m 

Interval  width  (DST)  fi m 

Out  values:  Number  of  unit  elements  (I) 

Partial  volume  for  the  various  unit  elements 
(DEL) 

Speed  of  descent  of  the  various  unit  elements 
(VZ)  m/s 

Calculations  Procedure 


The  droplet  size  distribution  in  the  cloud  element 
is  given  by  the  equation 


dV  -  56.V0  <  x5 
dX  5!  .D6 


exp 


l  L> 


/ 


in  which  X  is  the  droplet  diaiuete1'  in  /ira 
VQ  is  the  total  volume  of  substance  in  the  cloud 
element 

D  ■  .817  •  DMM 


The  value  of  the  total  volume  is  here  unimportant  and  is 
set  ■  1 


The  equation  is  solved,  giving: 


56 


exp 


51 


J 

o 


i  D 


X  }  dX  «  Vf 


•  exp  j -  1  X  f 
_L  '  D  ) 


5r  1 
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The  speed  of  descent  is  detarminee.  by 

VZ  =  k  •  X<C1+C2>  '  10 log  X  +  C3  •  (10log  X)2  ) 
k  =  3.0044  •  10  '5 

Cl  «  2,  C2  ■  C2  *  0  0  *  X  <  53  ft m 

,  Cl  .  2.01302,  C2  •  0.12055,  C3  -  -0.07422 
*  53  *  X  <  800  fm 


For  I  *  1,  2,  3  ...  .  calculate 

DEL(I)  ■  part  volume  for  the  interval  DST  •  I 

VZ(I)  ■  speed  of  descent  for  droplet  size  X  s 
DST  .  (2-1  -1) 


2 

Normally,  the  calculations  arc  broken  off  when  s.  part 
volumes  >.999,  whereupon  the  main  program  is  returned  to. 

If  this  condition  is  not  fulfilled,  the  calculations 
are  broken  off  when  the  number  of  unit  elements  •  100,  with 
the  end  line: 

"NUMBER  OF  UNIT  ELEMENTS  -  100.  REMAINING  VOLUME  XXX" 

with  return  to  the  main  program. 

For  a  normal  break  off  the  width  of  the  interval 
should  be  5s  10  4m, 
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*TA*r 


RETURN 


-  35  - 


Flow  Scheme,  Main  Program 

_/Legend7:  (1)  Table  of  4  as  a  function  of  wind 
strength  and  layering;  (2)  End  of  data  card?; 

(3)  Yes;  (4)  Read  in:;  (5)  Determine  row  No. 
ISN  A  table  for  layering;  (6)  Read  in:;  (7) 

No;  (8)  Determination  of  median  wind  strength; 
(9)  Determine  column  No.  J  in  A  table  for  wind 
strength  VM  (I);  (10)  Call  MOFHST  for  division 

un^  elements;  (11)  Read  in  and  write  out; 
(12)  Estimate  part  amount  substance  DW  (I)  for 
unit  element  I-NF;  (13)  Write  out  table  of 
division  into  unit  elements;  (14)  Write  out 
BLG  -  coverage  in  XM  (I),  YM  (I)  at  time  TS  (I). 


Flow  Scheme  MOFHST 

/Legend7:  (1)  Write  out:  "Number  of  Unit 
Elements  ■  100,  Remaining  volume  ■  LIV  /7sic7; 
(2)  No;  (3)  Yes.  X  “ 


-  END  - 
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